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A new cosmological model based on the de Sitter gravity is investigated by dynamical analysis 
and numerical discussions. Via some transformations, the evolution equations of this model can 
form an autonomous system with 8 physical critical points. Among these critical points there exist 
one positive attractor and one negative attractor. The positive attractor describes the asymptotic 
behavior of late-time universe, which indicates that the universe will enter the exponential expansion 
phase, finally. Some numerical calculations are also carried out, which convince us of this conclusion 
derived from the dynamical analysis. 
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Iii the last decades, some new cosmological observations such as SNela, CMBR and large scale structure all indicate 
that our universe is accelerating expanding and there exists a new mystical energy component in our universe, 
dubbed dark energy. [1] Many heuristic models have been proposed to explain the nature of this new component, 
which account for almost 74% of the energy density of our universe. Some of them have physical foundations, [2] 
some are just phenomcnological. [3] However, no one is flawless. All these models have their own problems, such as 
cosmological constant problem and fine tuning problem. Recently, a new kind of model, called torsion cosmology, has 
drawn researchers' attention, which is basically based on some new gravity theory, the gauge theory of gravity. Among 
these models, Poincare gauge theory is the one that has been investigated widely, which is inspired from the Einstein 
special relativity and the localization of Poincare symmetry. [4] Nester et al.[b] applied this new gravity theory to 
Oh' cosmology, and obtained a novel model which is likely to be a new explanation of the accelerating expansion. In that 
model, dynamic connection mimics the contribution of dark energy. Based on the work of Nester and his colleagues, 
some dynamics analysis and analytical discussion have been conducted in many papers, from which we can know the 
' fate of the universe more clearly. [6] 

Besides Poincare gauge theory, there is another classical gauge theory of gravity, de Sitter gravity, which can also 
be the alternative gravity theory for Einstein Gravity. This theory is derived from the de Sitter invariant special 
, relativity and the localization of de Sitter symmetry. [7] In the de Sitter gravity theory, Lorentz connection and 
' tetrad are combined to form a new connection, i.e. dS connection, which is valued in so(l,4), rather than Lorentz 
connection's so (1,3); and the gravitational action takes the form of Yang-Mills gauge theory. Like Poincare gauge 
theory, the spacetime also has a generic Riemann-Cartan structure, U±. From the variational principle one can obtain 
the gravitational field equation. Analogous to PG theory, de Sitter gravity has also been applied to the cosmology 
. i recently to explain the accelerating expansion. [8] 

In this Letter, we transform the cosmological evolution equations of de Sitter gravity model into a set of first-order 
dynamics equations, which form an autonomous system. Then, we give some dynamical analysis of this autonomous 
system and obtain all the critical points. We analyze the critical points' dynamical properties and stabilities, and 
. . , find out that there exists a late-time de Sitter attractor. It is concluded that the universe will expand exponentially 
k*" ■ in the end, as the ACMD predicted. 

i For a homogeneous, isotropic universe, the space-time is described by the FLRW metric: 

ds 2 = dt 2 - a 2 (t)[-^— + r 2 (d9 2 + sin 2 t?# 2 )], (1) 
1 — for 

and the isotropic and homogeneous torsion takes the form 

T° = 

T 1 = T + $° A tf 1 + T_ d 2 A tf 3 

T 2 = T+ i?° A d 2 - T_ d 1 A $ 3 (2) 
T 3 = T+ A tf 3 + T_ 1 A i? 2 , 
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where T + is the trace part of the torsion, \T a i dal while T_ is the traceless part of the torsion. They are all functions 
of time t, with + and — denoting the even and odd parities, respectively. Here #° = dt, i? 1 = a(t)dr, d 2 = 
a(t)rd6 and i? 3 = a(i) sin 6»d0. 

According to the field equations of dc Sitter gravity and Eqs.(l) and (2), one can easily obtain the new evolution 
equations of universe, [8] 
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Equations (3) and (4) are the 0-0 and 1-1 component of Einstein-like equations, respectively; and Eqs.(5)-(6) are two 
independent Yang-like equations, which are derived from the (r, 9, (f>) and (t, r, r) components. Here we also assume 
that the spin density is zero. 

Also, the energy momentum tenor is conserved by the virtue of Bianchi identities, leading to the continuity equation 



3d , . 
P = - — (P + P)- 



(7) 



Equation (7) can be derived from Eqs.(3)-(6), which means 4 of Eqs.(3)-(7) are independent. These four in- 
dependent equations with the EOS of matter content comprise a complete system of equations for five variables 
a(t), T_(t), T + (t), p(t) andp(i). By some algebra and differential calculations, we could simplify these five equations 
to 
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where H = a/a denotes the Hubble parameter. In order to make these equations dimensionless, we can rescale the 
variables and parameters as 



t -> t/l ; H -> / ff ; fc -)■ Zgfc; i? ->• #/Z ; 

T 4ttGZ 2 4ttGZ 2 

T±^l T±: p^ — — ^p; p^ — — ^p, 



(13) 
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where Iq = c/Hq is the Hubble radius. 

By some further calculations, we find that if the equation of state is constant, then these equations would turn out 
to form a six-dimensional autonomous system, which takes the forms 

2 / n \ 3 (i+™) 3 / T 2 \ 

H = -2H 2 + — 2 k \J-j +-l K P + 3HT + -Tl + —)+(l + 3w)p, (14) 



+6H(l + w)p, 
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where p is a dimcnsionless parameter denoting the current energy density. For such an autonomous system, i.e. 
Eqs.(14)-(19), we can use the qualitative method of ordinary equations with respect to the new set of variables, 
(if, P, Q, T + , T-, p). In order to find the critical points of this system, we should set the left-hand side of 
Eqs.(14)-(19) to zero and solve these algebra equations. Next, we just consider the matter dominant case with 



spatial flatness. We find that there are nine critical points {H c , P c , Q c 



T +C ,T- 



~, p c ) of this system, as shown in Table 1. 



Furthermore, we analyze the stability of these critical points by means of first-order perturbation. By the Taylor 
expansion, we could obtain the perturbation equation around the critical points, i.e. 

5x = Ax, A = ^-\ x=Xc , (20) 

where x means the six variables of this autonomous system and / denotes the six vector functions on the right-hand 
side of Eqs. (14)-(19). According to the dynamical analysis theory, we could classify these critical points by the 
coefficient matrix A's eigenvalue. The classification of these critical points are shown in Table 2. It is easy to find 
that there are only one positive attractor, i.e. point 1, whose eigenvalues all are negative, and only one negative 
attractor, i.e. point 2, whose eigenvalues all are positive. Positive attractors attract the solutions nearby it, while 
negative attractors repel them. The phase lines connecting a positive attractor and a negative attractor are called 
the heteroclinic lines as shown in Fig. 1. 

Critical points are actually some exact solutions of a dynamic system, especially the positive attractors which are 
often the extreme points of the orbits in phase space. Therefore, they would describe the asymptotic behaviors of 
solutions. In this model, the positive attractor point 1 would show us the late time evolution of our universe. It is 
indicated that all quantities tend to zero, except the Hubble constant, which will remain at a fixed value, and therefore 
the whole universe will approach the exponential expansion, just like the ACDM model. 

Hoo = - a = i =► a(t) cx cxp (*) (21) 



As we know, the critical points can only describe the local properties. If we want to know the global properties, 
we have to resort to the numerical calculations. We solve Eqs.(14)-(19) numerically, and show some generic solutions 
in Fig. 2. From these numerical results, we could see easily that the late-time positive attractor covers a wide range 
of initial values, and therefore alleviate the fine-tuning problem. Also for its insensitivity to the initial conditions, 
this autonomous system has no chaotic features. From this perspective, the cosmology based on de Sitter gravity is 
quite different from the one based on the PG theory, which suggests that the expansion will asymptotically come to 
a halt. [6] However, this discrepancy is only due to the existence of the de Sitter radius R. If we set R — > oo, the de 
Sitter gravity would degenerate to the PG theory, and have the same conclusion. 

In summary, we have studied the torsion cosmology based on de Sitter gravity. According to Ref. [8] , we have rewrit- 
ten the evolution equations as a set of first-order dimensionless equations which form a six-dimensional autonomous 
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Critical points Eigenvalues 

(1) (i, 0,0, 0,0,0) 

(2) (-i, 0,0, 0,0,0) 

(3) (-^, 0,0,-|, 0,0) 

(4) (^,0,0,|,0,0) 

(5) (-^,0,0,^,0,0) 

(6) (^ 0,0,-^,0,0) 

(7) (0,0,0,-^,0,^) 



(8 ) (0,0,0,^,0,^) 



(9) (0,0,0,0,0,^1) 



TABLE I: The critical points and their corresponding eigenvalues. The point 9 is not physical, for its negative energy density. 



Critical points 


Property 


Stability 


(1) 


Positive-attractor 


Stable 


(2) 


Negative-attractor Unstable 


(3) 


Saddle 


Unstable 


(4) 


Saddle 


Unstable 


(5) 


Saddle 


Unstable 


(6) 


Saddle 


Unstable 


(7) 


Spiral-saddle 


Unstable 


(8) 


Spiral-saddle 


Unstable 



TABLE II: The stability properties of critical points. 

system. We find out that among all the eight physical critical points, there are one positive attractor and one negative 
attractor. The positive attractor implies that the universe will expand exponentially in the end and all other physical 
quantities will turn out to vanish. Also we present some numerical analysis of this model, and find out that the 
late-time evolution is not sensitive to the initial values and parameter and, for a large range of parameter choice, the 
dynamical system would approach to the positive attractor. Therefore, in this sense, the de Sitter gravity model looks 
more like the ACDM model, [2] rather than the PG theory. [6] 

If we want to know deeper on whether this model can explain the accelerating expansion, we have to settle the initial 
values and parameter choice. It requires us to do some further analytical and numerical calculations and examine 
them with the current observations, such as SNela and CMB etc. These issues will be considered in the upcoming 
papers. 
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FIG. 1: The (H,T+,p) section of the phase diagram with R — 4/3. The heteroclinic orbits connect the critical points 1 and 2. 
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FIG. 2: Evolution of Hubble constant H with respect to some initial values and parameter choice (R, Hq, Po,Qo,7+o, 

T_o, po). According to the transformations (13), the unit of time here is the Hubble Time, (a) R is fixed and T± is changed. 

(b) R is changed. 
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